T = Phase Diagram of the Double-Exchange Model 
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We present the T — phase diagram of the double-exchange model (ferromagnetic Kondo lattice 
model) for all values of the carrier concentration n and Hund's couplng J, within dynamical mean 
field theory. We find that depending on the values of n and J, the ground state is either a ferromag- 
net, a commensurate antiferromagnet or some other incommensurate phase with intermediate wave 
vectors. The antiferromagnetic phase is separated by first order phase boundaries and wide regimes 
of phase separation. The transition from the ferromagnetic phase to an incommensurate phase is 
second order. 



I. INTRODUCTION 

In this Brief Report we present the T = CUmagnetic 
phase diagram of the double exchange modeltl believed 
relevant tO|-,the physics of 'colossal' magnetoresistance 
manganitesa. Our results may also he-of interest in the 
context of magnetic semiconductors!]'!! We obtain the 
different magnetic phases, the boundaries between them, 
and the regimes of phase separation using the dynamical 
mean field approximation (DMFT)ETlI We consider only 
the double-exchange interaction, and study only mag- 
netic ordering, not charge or orbital ordering. This pa- 
per is a companion to a previous papercl which reported 
results mainly in regimes believed directly relevant to ex- 
periment on manganites. Here we present the complete 
T = phase diagram of the model, for the convenience 
of other workers, and to provide a basis for quantitative 
comparison of the dynamical mean field method to other 
methods, i.e the phase diagram given in this paper is 
not just limited to CMR physics in manganites, but is of 
much more general validity covering other systems like 
Mn doped GaAsQ. 

The double-exchange-only (DE) model involves elec- 
trons moving in a band structure defined by a hopping 
matrix tfj and a chemical potential /z, and connected by 
a Hund's coupling (J > 0) to core spins S. We denote 
the operator creating an electron in orbital a with spin 
a on site i by d\ aa and define the double-exchange only 
Hamiltonian Hoe b; 



h de =- t % d i a d ]ba -^ d L d - 

<ij>a ia. 

— J ^ S c .i • d ia a aj 3dip . 



(i) 



where < ij > implies sum over nearest neighbors. 

Note that there are only two independent dimension- 
less parameters defining Hoe, since the hopping ampli- 



tude can be used to define the energy unit. Our interest 
in this paper is the T = magnetic phase diagram of 
H d e as a function of density n and coupling J. To com- 
pute this we use the dynamical mean field approximation 
which becomes exact in the limit of spatial dimensional- 
ity d — > oo and is believed to be reliable in d = 3 [Ref. 
|9| . The complete T — phase diagram of this important 
model has not to ouc, [knowledge been presented in any 
limit. Other workersSllj have considered two phases: a 
commensurate antiferromagnet (AF) and a homogeneous 
ferromagnet (FM) . In a previous papera we demonstrated 
the existence of a free energy minimum with spin or- 
der which is neither AF nor FM; this possible phase has 
magnetic correlations incommensurate with the underly- 
ing lattice and was denoted IC. However, in our previous 
work the precise spin correlations were not established, 
nor was the issue of phase separation addressed. The 
importance of phase separation p,Ae CMR context has 
been stressed by other workersU'EI. Here, we resolve 
these issues. As results obtained by other methods be- 
come more available, comparison to our results will allow 
the accuracy of DMFT in d — 3 to be assessed. 

For our specific computations we take the core spins 
to be classical. This is reasonable because we consider 
T = energetics in high spatial dimension, and in cases of 
physical relevance the core spin is not small (S = 3/2 for 
the manganites). We proceed by calculating the energies 
of the different phases. We identify the ground state as 
the phase of lowest energy and use an explicit Maxwell 
construction to determine regions of phase separation. 

We now discuss the possible phases—which may occur. 
In the d — > oo limit, Muller-Hartmannlij pointed out that 
the susceptibility depends on the one-dimensional param- 
eter X(k) = 2^2i cos (ki)'-> similar considerations imply 
the ordering pattern and the energy are also functions of 
X. Because — 1 < X < 1, we write X — cos((9); 6 = 
corresponds to ferromagnetic (FM) phase and 8 = it to 
antiferromagnetic (AFM) phase considered previously^. 
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< 9 < 7r coresponds to a phase with incommensurate 
(IC) spin correlations. The specific IC phase treated in 
our previous workQ corresponds to — n/2. 

The DMFT solution of Eq. @ may be written in terms of 
a spin-dependent mean field function a(cj), which obeys 
an equation that depends on the density of states of the 
underlying lattice. We study the semicircular density 
of states N(e) — y/At 2 — e 2 /(2-7rf 2 ), corresponding to the 
d — > oo limit of the Bethe lattice with nearest neighbor 
hopping, so that a^w) = lo + fi - J2s t s Gl i+s,i+s t s witn 
tg the hopping. To solve the equations, we adopt a space 
dependent spin basis, quantized parallel to the core spin 
on site i. The conventions just described then imply that 
a has two independent components aj , a j. given by 

, , , cos 2 ((9/2) 2 sin 2 (6>/2) 

. . 2 cos 2 (6>/2) 2 sin 2 (672) 

aj,(wj — J a,\[io) + J 

These equations contain as special cases the equations 
previously studied for the FM, AFM and 9 = n/2 phases. 
They may be obtained by minimizing with respect to 
a T ,ax, the energy 



E = TJ2 



[a(u n ) - (u„ + n)] 2 b(uj n f 



t 2 t 2 cos9 
■ln[oK) a -(6(w n ) + J) 2 ] (3) 



where a = (o-f +ajJ/2 and & = (o-f — oj.)/2. This equation 
may be transformed into the familiar expressions 

//•OO 7 
rfe fc ^V(e fe ) y — /(w)wTr[A(e fc ,w)] 

/OO J 
-/( W ) W Tr[ImG ioc ( W )] (4) 
-oo T 

through an integration by parts and use of the mean field 
equations (recalling that E is extremized by a-f, aj,). Here 
N{ek) is the "band" density of states determined from the 
hopping tfj, A(e k ,uj) = — ^Im.G(ck, w) and the electron 

Green function G = (w+,u— e^— S(w)) -1 . S = a— G ; ~J, is 
the self energy, which is local (k-independent) but depen- 
dent on spin in the dynamical mean field approximation. 
Replacing the form of G; oc in Eq. |[ we get 



dto ( 2a(uj) 



(5) 



We have solved the mean field equations (Eq. ||) and 
computed the energies for a wide range of 9, n. The re- 
sulting phase diagram is shown in Fig. [j]. As one can see 
from Fig. [I], a homogeneous ferromagnetic phase is sup- 
ported over a wide range of fillings and couplings. There 
is a small island of IC phase characterized by a 9 which 
varies with position inside the IC phase. The AFM phase 



exists only at n = 1. The FM-IC phase boundary is ap- 
prently second order; all of the other phase boundaries 
are (within our numerical accuracy) first-order and are 
accompanied by regions of phase separation. At suffi- 
ciently small J, phase separation preempts the IC phase 
entirely even though this phase minimizes the energy for 
0.47 < n < 0.7. 

Fig. H shows energy as a function of n for the FM 
and AF states, as well as 9 = ir/2 and 3tt/A. We see 
immediately from these curves that at n = 1 the AF 
phase is the lowest; at intermediate n there is a sequence 
of IC phases characterized by a varying angle and at small 
n the FM phase is stable. Between the n = 1 AF phase 
and the intermediate n IC phase there is a region of phase 
separation, and the 9 value at the PS-IC boundary is 
a non-universal J-dependent number. Fig. ^ shows the 
behavior in the vicinity of the FM/IC phase boundary. 
We see that the transition is apprently second-order, with 
9 evolving smoothly to at the FM/IC phase boundary. 

Fig. [| shows the behavior for J = 0.6£ and n close 
to 1. A Maxwell construction shows that the 9 = 7it/8 
curve lies higher than the line connecting the AF n = 1 
point and the 9 — 3tt/A, n sa 0.78 point. The energy 
differences involved are seen to be tiny, implying large 
uncertainties in the location of the phase boundary and 
rather complicated behavior in physical systems. 

A generally noteworthy feature of Fig. |^ is the ex- 
tremely small energy differences between different phases: 
typically the energy difference per site AE between the 
lowest and next lowest phase is only aboutJ0.02i. For 
manganites, a reasonable value of t is 0.65eVcl leading to 
AE/site ~ 100K. This suggests an unusuaLsensitivity 
to disorder, as already noted by Moreo et al.li3. 

We now use analytic arguments to analyze more care- 
fully the behavior in the small J limit, where the energy 
may be written as 



E = -\ X {0)J 2 



with 



X (9) = -2 J ^/Hlm 
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ao(w) 2 — t 2 cos 6* 



(6) 



(7) 



and a = g [ w + A* — v( w + m) 2 ~ 4>i 2 ] 

Note that the behavior as 9 — > ir and n — > 1 is com- 
plicated because the energy of the n = 1,6 = w AFM 
phase ~ J 2 ln(l/J). From the explicit energy expression 
we find (neglecting 0(J 4 )) 

9 -§f = {[E ( Sl n~^/2); -4/ J 2 ) - E (vr/2; —4/ J 2 )] - 



4J 



[F {sin-^/2); -A/ J 2 ) - F (tt/2; -4/ J 2 )] (8) 



where F((f>; m), E(<fi; m) are elliptic integrals of the First 
and Second Kind. The corresponding energies are plotted 
in Fig. H for J = 0.05. Note that for n very near 1 
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(fx < J), \af is linear in 1 — n. Thus at n ^ 1 the 
AF phase is always unstable to phase separation. The 
Maxwell construction also shows that in the small J limit 
the IC phase is preempted entirely, and the ferromagnetic 
phase is only stable for n < 0.15. 

We now briefly discuss the physical meaning of the IC 
phase. From Eq. ^], the lowest energy phase is the one 
with the largest x(6>). For typical band structures even 
in d = 3 x takes one value at q = and another value 
at q = Q = (7r,7r,7r); the variation between these lim- 
its is smooth for momenta along the zone diagonal, but 
over much of the Brillouin zone \ is quite weakly momen- 
tum dependent. In the d — ► oo limit this weak momen- 
tum dependence disappears: at any q ^ A(7r,7r, ...) 
takes the g-independent value \ioc = ^2 q x{o) an d along 
the diagonal. For d large but finite, x varies rapidly 
with ^crossing over from its q = X(tt,tt,...) values to a 
very weakly g-dependent value in a range Sq « 0(l/y/d). 
Thus the hallmark of the phases which are characterized 
by a small cos 8 is a large range of spin arrangements very 
nearby in energy. We suspect that this regime would ap- 
pear experimentally as a spin glass in the presence of 
disorder 

In conclusion, we have determined the T = phase 
diagram of the double-exchange Hamiltonian including 
regions of phase separation, using analytic arguments at 
small and large J and the dynamical mean field method- 
We find that the incommensurate phase noted in Ref.El 
has a rather small region of existence, being mostly over- 
whelmed by phase separation, and have found the differ- 
ences in energy between the different phases to be very 
small. 

We have learned of related work by L. Yin and J. Ho 
(APS March Meeting 2000 Bulletin; L. Yin and J. Ho, 
private communication) , who used a different mean-field 
technique in d = 3 to obtain a phase diagram including 
spiral states, with regions of phase separation qualita- 
tively similar to what is found here. Our results are also 
similar to the phase diagram of Alonso et al£3. 
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FIG. 1. The phase diagram of the double exchange model, 
as deduced from the analytic calculations of the free energy 
at all J. 
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FIG. 2. Energy per orbital for J/t = 0.6. There is a 
continuity of phases from 9 = to 8 — 3n/4, followed by 
phase separation between 9 = 3n/4 and 9 = n 
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FIG. 3. The region between 0.2 < n < 0.6 has been ampli- 
fied to illustrate that the transitions between 9 — to 9 = n /2 
are second order. The second order transitions extend to the 
e = 3tt/4 state 




-l.l 1 • 1 • ' 

0.75 0.85 0.95 



n 

FIG. 4. The transition from the 9 = 37r/4 state to the 
AFM state is first order, and the system phase separates. 




FIG. 5. Small J energies of the states with different 9 
computed from Eqs. and [] at J = 0.05t. The wide regime 
of phase separation between 9 = 0, n — 0.15 and 9 = n, n = 1 
phases is evident. 
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